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Retransmission-based failure recovery represents a primary approach in existing com- 
■ munication networks that guarantees data delivery in the presence of channel failures. 

Recent work has shown that, when data sizes have infinite support, retransmissions can 
cause long (-tailed) delays even if all traffic and network characteristics are light-tailed. In 
this paper we investigate the practically important case of bounded data units < Lf, < 5 
under the condition that the hazard functions of the distributions of data sizes and channel 
^ \ statistics are proportional. To this end, we provide an explicit and uniform characteriza- 

tion of the entire body of the retransmission distribution P[iVb > n] in both n and b. Our 
main discovery is that this distribution can be represented as the product of a power law 
and Gamma distribution. This rigorous approximation clearly demonstrates the coupling 
of a power law distribution, dominating the main body, and the Gamma distribution, de- 
04 ■ termining the exponential tail. Our results are validated via simulation experiments and 

\l I can be useful for designing retransmission-based systems with the required performance 

characteristics. From a broader perspective, this study applies to any other system, e.g., 
\ computing, where restart mechanisms are employed after a job processing failure. 



1 Introduction 



Failure recovery mechanisms are employed in almost all engineering systems since complex 
d \ systems of any kind are often prone to failures. One of the most straightforward and widely 

used failure recovery mechanism is to simply restart the system and all of the interrupted 
jobs from the beginning after a failure occurs. It was first recognized in [6l [18] that such 
mechanisms may result in long-tailed (power law) delays even if the job sizes and failure rates 
are exponential. In [12], it was noted that the same mechanism is at the core of modern 
communication networks where retransmissions are used on all protocol layers to guarantee 
data delivery in the presence of channel failures. Furthermore, [12] shows that the power 
law number of retransmissions and delay occur whenever the hazard functions of the data 
and failure distributions are proportional. Hence, power laws may arise even if the data and 
channel failure distributions are both Gaussian. In particular, retransmission phenomena can 
lead to zero throughput and system instabilities, and therefore need to be carefully considered 
for the design of fault tolerant systems. 



This work is supported by the National Science Foundation grant number 0915784. 
Extended abstract of this paper with preliminary results has appeared earlier in [8] . 
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More specifically, in communication networks, retransmissions represent the basic building 
blocks for failure recovery in all network protocols that guarantee data delivery in the presence 
of channel failures. These types of mechanisms have been employed on all networking layers, 
including, for example. Automatic Repeat reQuest (ARQ) protocol (e.g., see Section 2.4 of [4]) 
in the data link layer where a packet is resent automatically in case of an error; contention based 
ALOHA type protocols in the medium access control (MAC) layer that use random backoff 
and retransmission mechanism to recover data from collisions; end-to-end acknowledgement for 
multi-hop transmissions in the transport layer; HTTP downloading scheme in the application 
layer, etc. It has been shown that several well-known retransmission based protocols in different 
layers of networking architecture can lead to power law delays, e.g., ALOHA type protocols 
in MAC layer [13^ ITO] and end-to-end acknowledgements in transport layer [11^ [9] as well as 
in other layers [12]. For other (non-retransmission) mechanisms that can give rise to heavy 
tails see |15j and the references therein. In particular, the proportional growth/multiplicative 
models can result in heavy tails [151 E] ■ 

Traditionally, retransmissions were thought to follow light-tailed distributions (with rapidly 
decaying tails), namely geometric, which requires the further assumption of independence 
between data (packet) sizes and transmission error probability. However, these two are often 
highly correlated in most communication systems, meaning that longer data units have higher 
probability of error, thus violating the independence assumption. Recent work dH [9] 

has shown that, when the data size distribution has infinite support, all retransmission-based 
protocols could cause heavy-tailed behavior and possibly result in zero throughput, regardless 
of how light-tailed the distributions of data sizes and channel failures are. Nevertheless, in 
reality, packet sizes are upper bounded by the maximum transmission unit. For example, 
WaveLAN's maximum transfer unit is 1500 bytes. This fact motivates us to investigate the 
transmission of bounded data and approximate uniformly the entire body of the resulting 
retransmission distribution as it transits from the power law to the exponential tail. 

We use the following generic channel with failures [T2] to model the preceding situations. 
This model was first introduced in [6j in a different application context. The channel dynamics 
is described by the i.i.d. channel availability process {A, Aj}j>i, where the channel is contin- 
uously available during periods {Ai} and fails between these periods. In each period of time 
that the channel becomes available, say Ai, we attempt to transmit the data unit of random 
size Lf,. We focus on the situation when the data size has finite support on interval [0,6]. If 
Lb < Ai, we say that the transmission is successful; otherwise, we wait for the next period ^j+i 
when the channel is available and attempt to retransmit the data from the beginning. It was 
first recognized in [6] that this model results in power law distributions when the distributions 
of L = Loo and A have a matrix exponential representation, and this result was rigorously 
proved and further generalized in [121 El E]. A related study when L = ^ is a constant and 
failure/arrival rates are time-dependent Poisson can be found in [3]. 

It was discovered in p2] that bounded data units result in truncated power law distribu- 
tions for the number of retransmissions, see Example 3 in [12j . see also Example 2 in |13j . Such 
distributions are characterized by a power law main body and an exponentially bounded tail. 
However, the exponential behavior appears only for very small probabilities, often meaning 
that the number of retransmissions of interest may fall inside the region where the distribution 
behaves as a power law. It was argued in Example 3 of that the power law region will 
grow faster than exponential if the distributions of A and L), are lighter than exponential. The 
retransmissions of bounded documents were further studied in [19], where partial approxima- 
tions of the distribution of the number of retransmissions on the logarithmic and exact scales 
were provided in Theorems 1 and 3 of [19j, respectively. In this paper, we present a uniform 
characterization of the entire body of such a distribution, both on the logarithmic as well as 
the exact scale. 
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Specifically, let Nfj represent the number of retransmissions (until successful transmission) 
of a bounded random data unit of size Lh G [0, b] on the previously described channel. In order 
to study the uniform approximation in both n and b we construct a family of variables Lh, 
such that ¥[Lb < x] = P[L < x]/F[L < 6], for < x < 6 when L = Loo is fixed. This scaling of 
Lfe was also used in [19j. For the logarithmic scale, our result stated in Theorem 12.21 provides 
a uniform characterization of the entire body of logP[A'^f, > n], i.e., informally 



for all n and b sufficiently large when the hazard functions of L and A are linearly related 
as logP[L > x] ~ alogP[74 > x]; see Theorem 12.21 for the precise assumptions. Note that 
the first term in the preceding approximation corresponds to the power law part n~" of the 
distribution, while the second part describes the exponential (geometric ¥[A < 6]") tail. Hence, 
it may be natural to define the transition point Uf, from the power law to the exponential tail 
as a solution to nblogF[A < 6] « alognj,. 

In addition, under more restrictive assumptions, we discover a new exact asymptotic for- 
mula for the retransmission distribution that works uniformly for all large n, b. Surprisingly, 
the approximation admits an explicit form (see Theorems 12.31 and 12. 4p 



where x Ay = min(x, y) and £{■) is a slowly varying function; note that the preceding integral 
is the incomplete Gamma function T{x, a). 

Clearly, when — nlogP[^ < i 0, the preceding approximation converges to a true power 
law T{a + l)/{i{n)n"). And, when — nlog(P[A < b]) t oo, approximation (jl.ip . by the property 
T[x,a) ~ e~^x"~^ as x — )• oo, has a geometric leading term ¥[A < 6]". Interestingly, for the 
special case when a is an integer and l{x) = 1, one can compute the exact expression for 
P[A'^E, > n], see Proposition 12.21 Furthermore, our results show that the length of the power 
law region increases as the corresponding distributions of L and A assume lighter tails. All of 
the preceding results are validated via simulation experiments in Section [3l It is worth noting 
that our asymptotic approximations are in excellent agreement with the simulations. 

This uniform approximation allows for a characterization of the entire body of the distribu- 
tion ¥[Nb > n], so that one can explicitly estimate the region where the power law phenomenon 
arises. Introducing the relationship between n and P[^ > b] also provides an assessment method 
of efficiency and is important for diminishing the power law effects in order to achieve high 
throughput. Basically, when the power law region is significant, it could lead to nearly zero 
throughput (a < 1), implying that the system parameters should be more carefully adjusted 
in order to meet the new requirements. On the contrary, if the exponential tail dominates, the 
system performance is more desirable. Our analytical work could be applicable in network pro- 
tocol design, possibly including packet fragmentation techniques [HI |T6j and failure-recovery 
mechanisms. In particular, since we approximate the entire body of the distribution, we could 
estimate the mean value of A'^;,! and thus, via Wald's identity, the mean delay and throughput. 

Also, from an engineering perspective, our results further suggest that careful re-examination 
and possible redesign of retransmission based protocols in communication networks might be 
necessary. Specifically, current engineering trends towards infrastructure- less, error-prone wire- 
less technology encourage the study of highly variable systems with frequent failures. Hence, 
our results could be of potential use in improving the design of future complex and failure-prone 
systems. 

The rest of the paper is organized as follows. After a detailed description of the channel 
model in the next Subsection 11.11 we present our main results in Section [2j Then, Section [3] 



logP[A'b > n] 



a log n + n log P[^ < b] 




(1.1) 
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contains simulation examples that verify our theoretical work. A number of technical proofs 
are postponed to Section [H 



1.1 Description of the Channel 



In this section, we formally describe our model and provide necessary definitions and notation. 
Consider transmitting a generic data unit of random size Li, over a channel with failures. 
Without loss of generality, we assume that the channel is of unit capacity. As stated in the 
introduction, the channel dynamics is modeled by the channel availability process {A, Ai}i>i, 
where the channel is continuously available during time periods {Ai} whereas it fails between 
such periods. In each period of time that the channel becomes available, say Ai, we attempt to 
transmit the data unit and, if Ai > Lf,, we say that the transmission was successful; otherwise, 
we wait for the next period Aj+i when the channel is available and attempt to retransmit the 
data from the beginning. A sketch of the model depicting the system is drawn in Figure [H 



Failure-prone 
channel 

{An} 




resend no 



Figure 1: Packets sent over a channel with failures 

We are interested in computing the number of attempts Nf, (retransmissions) that it takes 
until Lf) is successfully transmitted, which is formally defined as follows. 

Definition 1.1 The total number of retransmissions for a generic data unit of length Lfy is 
defined as 

Nb ^ inf{n : A^ > U} 

We denote the complementary cumulative distribution functions for A and L, respectively, as 

G{x) = ¥[A > x] 

and 

F{x) =¥[L > x], 

where L is a generic random variable that is used to define the distribution of Lf,. 

Throughout the paper we assume that L and A are continuous (equivalently, F{x) and 
G{x) are absolutely continuous) and have infinite support, i.e., G{x) > and F{x) > for all 
X > 0. Then, the distribution of Lf, is defined as 



[Lb < x] 



¥[L < x] 



0<x<b. 



X2) 



F[L < b] ' 

To avoid trivialities, we assume that b is large enough such that P[L < b] > 0. 

In this paper we use the following standard notations. For any two real functions a{t) and 
b(t) and fixed to S M|J{oo}, we use a{t) ~ b{t) as t — )• to to denote limt^to CL{t) /b{t) = 1. 
Similarly, we say that a{t) > b{t) as t — )• to if hminft^to [a(t)/6(t)] > 1; a(t) < b{t) has a 
complementary definition. 
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2 Main Results 



In this section, we present our main results. Under mild conditions, we first prove a gen- 
eral upper bound for the distribution of Ni, on the logarithmic scale in Proposition 12. li In 
Theorem 12.21 present our first main result, which under more stringent assumptions, char- 
acterizes the entire body of the distribution on the logarithmic scale uniformly for all large n 
and b, i.e., informally we show that 

log F[Nb > n] -alogn + n log¥[A < b] , 

as previously mentioned in the introduction. Roughly speaking, when — logP[A < b] = 
o(logn/n), ¥[Nh > n] is a power law of index a. Our results on the exact asymptotics are given 
in the next Subsection l2.1l in Theorems 12.31 and 12. 4t the results are stated in two different the- 
orems since Theorem 12.41 requires slightly stronger assumptions. The uniform approximation 
implied by these two theorems is presented in (|2.6p , or previously in (|l.ip . 

Recall that the distribution of Lf, has finite support on [0, 6], given by (|1.2p . First, we prove 
the following general upper bound. 

Proposition 2.1 Assume that 

logP[L > x] < alogP[A > x] as X oo 

and let bo be such that P[L < b^] > 0, < b^] > 0, then for any e > 0, there exists no, such 
that, for all n > nQ,b > bo, 

logP[iVfe > n] < (1 - e) [nlogP[^ < b] - a log n] . 

Remark 1 Note that this result can be restated as 

¥[Nb >n]< F[A < 6]"(i-^)n~"(i-^), 

for n, b sufficiently large. Hence, the distribution ¥[Nb > n] is bounded by the product of a 
power law and a geometric term. 

Proof: By assumption, there exists < e < 1 such that for all x > x^ > 6o > 0, 

F(x) < G(x)"(^-^). (2.1) 

Next, it is easy to see that F[Nb > n\Li,] = (1 — G{Lh))^, and thus, 

F[Nb >n]= E[l - G(Lb)]" 

= E[l - G(Lb)]"(i~^+^) 

< (1 - G(6))"(i-^) [E[l - G{Lb)rn{Lt < xe) + E[l - G(Lfe)]"^l(Lfe > x,)] 

dF{x) 



< (1 - G(6))"(i-^) 

< (1 - G(6))"(i-^) 



(1-G(x,)r+ / (l-G(x))' 



Fib) 







rj::+ / 1 -F(x)<^a-) 



dF{x) 



where rj^^ = 1 — G(xe), and the last inequality follows from (j2.ip : in case x^ > b > bo, 
the integral in the second inequality is zero and the last inequality trivially holds. Now, by 
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extending the preceding integral to oo, we obtain 



(1 - 



vZF{b)+ / l-F(x)^ dFix) 



_ I \ne 

<F(6)+E(1-F(L)^ 



where we use the elementary inequality 1 — x < e ^, x > 0, and thus 



mb>n] < (1 - 



<F(6)+Ee- 



_[/a(l-e) , 



by -F(-L) = U, where U is uniformly distributed on [0, 1] by Proposition 2.1 in Chapter 10 of 
[I7j : see also the Appendix. Hence, 



W>n]<-^(1-G(6)r(i-) 
(1 - 







(ne)"(i-^) 
q(1-6) 



(1 - (5(6)) 



n(l-e) 



<^(^)+(^|^r(a(l-e)) 



which follows from the definition of the Gamma function T{a) = Jq°° e ^dt. Therefore, 

P[iV6 > n] < (1 - G(6))"(i-^) 



< + ^.l^J-. r(a(l - e)) 



F(6)(ne)"(i-^)' 



< (1 - ^(b))"^!-^) 
= (1 - G'(6))"(^-^) 



C + ^^W4^^7r^r(a(l - e)) 



since b > bo, whereas, in the last inequality, we set H, = a{l - e)e-"(i-^)r(a(l - e))/F(6o). 
Now, we can choose tiq, such that for any e > and for all n > no, r/"^ < eH^n^°'^^^^\ so that 



F[Nb > n] < (1 - G(6))"(i-^) 



e — rr^ + 
n 



a{l-e) 



and by taking the logarithm in the preceding expression, we obtain 

logP[iVfc > n] < log {H,{1 + e)) + n(l - e) log(l - (5(6)) - a(l - e) logn 
= log(i/,(l + e)) + (l-e) [nlog(l - 0(6)) -Q logn] . 
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Next, since — nlog(l — G{b)) > and alogn > 0, n > 1, 



logP[iVfe >n] ^ log(ge(l + e)) _ (1 _ e) 



— nlog(l — G{b)) + alogn — nlog(l — G{b)) + alogn 

^ log(ffe(l + 6)) _ (1 _ ^ 

~ a log n 

and alogn being increasing in n, we can choose no such that for any n > no, 

log(g,(l + e)) 
alogn 

Thus, 

logPfiVfc > nl , , , , 

^ ^ ^ < -(1 - 2e), (2.2) 



< e. 



— n log(l — G{b)) + a log n 
which completes the proof by replacing e with e/2. □ 

Next, we determine the region where the power law asymptotics holds on the logarithmic 
scale. 

Theorem 2.1 If 

log F[L > x] ~ a log P[yl > x] as x ^ oo, (2.3) 

a > 0, then, for any e > 0, there exists positive uq, such that for all n > uq, for which 
n^+<^F[A > 6] < 1, we have 

-^-^nN,>-] - 1 < (2.4) 
alogn 

Note that this result appeared in Theorem 1 of [19]; for reasons of completeness, we present 
a simple proof in Section HI 

This result holds in the region n^^'' < 0{1/G{b)). Also, note that one can easily characterize 
the logarithmic asymptotics of the very end of the exponential tail of WlNf^ > n] for small b 
and large n. In particular, for fixed b, it can be shown that logP[A'^^ > n] ~ nlog(l — G{b)) as 
n — >• oo, see Theorem 1 in [19j. However, our objective in this paper is to determine the entire 
body of the distribution of 'S'[Ni, > n] uniformly in n and b. 

Next, we extend Theorem 12.11 to the entire region n > no, 6 > b^, which includes the 
geometric term ¥[A < 6]". For this theorem, we need slightly more restrictive assumptions. 
The reason why this is the case is that ¥[Ni) > n] behaves like a power law in the region where 
n = o(log n/G'(6)), while for n >> logn/G'(6), it follows essentially a geometric distribution; 
see Theorem 12.21 below. Hence, more restrictive assumptions are required since the geometric 
distribution is much more sensitive to the changes in its parameters (informally, ((l + e)x)~" ~ 
x-" but e-(i+^)^ 96 e-^). 

Definition 2.1 A function £{x) is slowly varying if i{x)/£{Xx) — )• 1 as x ^ 00 for any fixed 
A > 0. 

We also assume that functions i{x) are positive and bounded on finite intervals. 

Theorem 2.2 IfF[L > x] = i{¥[A > x]-'^)F[A > xf, for a > 0, i{x) slowly varying, then 
for any e > 0, there exist no, ^O; such that for all n > nQ,b > bo, 



■logP[iVfe>n] ^ 



-n log F[A < b] + a log n 



< e. (2.5) 
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Remark 2 Note that the statement of this theorem can be formulated in an equivalent form 



■logP[7Vb > n] 



n¥[A > 6] + a log n 
since —nlog¥[A < 6] ~ nF[A > 6] as 6 — ;> oo. 



< e, 



Remark 3 This theorem extends Theorem 1 in [19]. In particular, it proves the result uni- 
formly in n and b, while Theorem 1 in [T9j characterized the initial power law part of the 
distribution (n < G{b)~'^,0 < 77 < 1) and the very end with exponential tail (fixed b, 00). 

Proof: First, we prove the result in the region nG{b) > (51ogn, for any fixed 5 > 0, where 
the distribution of P[A'^[, > n] approaches a geometric tail. 

Note that, for all b > bo, the upper bound follows from Proposition 12.11 For the lower 
bound, observe that 



[Nb >n]= E[l - G(Lfe)]" 



Fix)-. 



> / 1 

1X0 \ 



\ dFjx) 

^^"(^(x)-!) J F{b) 

1 \ n 

F{x)^ 



dF{x) 



by F{b) < 1. Then, by the continuity of G{x), we can choose xq such that G{xo) = AeG'(6) 
where = A/(l — e), and A > 2 is such that G'(xo) < 1. Now, recall the slowly varying 
property of £{x) in the region {G{b)~^ , X~^G{b)~^} and thus, for bo large enough, b > bo, 



mb>n]> f (l-(l-e) ^("^^ 
Jxo \ 



E 



l-(l-6) 



F{Ly 



dF{x) 



£y^{G{b) 



— l{F{b)<F{L)<F{xo)) 



> 



{l-e)F{f))i/Q/£l/"(G{f))-l) 



(1 



which follows from F{L) = U, where U is uniformly distributed in [0, 1], and the change of 
variables z = {1 — e)U^^°' /£^^°'{G{b)^^). Thus, using our main assumption that F{bY/°' = 
£^/'^{G{b)-^)G{b), and F{xo)^/'^ = £^/'^{G{xo)-^)G{xo) = £^/'^{X:^G{by^)X,G{b) > (1 - 
e)AeG(6)^i/"(G(6)-i) = \G{b)£^/'^{G{b)-^), we obtain 

a£(GibV^) 

mb >n\> / _ (1 - zf z^-'dz. 



1-e 



(l-e)G(b) 



Now, if a > 1, is monotonically increasing and thus 



> aG{b) 



1 - e)° 
1 - zY+^ 



n + l 



J{l~e)G{b) 
(l-e)G{b) 

\(l-e)G(b) 
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aG{b) 



a(l+e)-l 



n + 1 



(1 - (1 - e)G(6))"+^ -(1 - A(l - e)G(6))"+^] 

n+l' 



n+1 



1 



l-A(l-6)G(b) Y 
1 - (1 - e)G{b) J 



where, for the second inequahty, we use the standard property of slowly varying functions (see 
[5]) that i{x) > x'""" for large x, and thus i{G{b)~^) > (5(6)'^", whereas the last inequality is 
implied by a > 1 and n + 1 < 2n for n > 1. 

Next, using the inequalities e~^^ < (1 — a^) < for x small enough, we have 



l-A(l-e)G(6) \ 

i-(i-6)G(6) ; 



n+l 



A(l-e)(n+l)G{6) 
< < e-(A-2){l-e)nG(fe) 



g-2(l-e)(n+l)G(6) 

< e~*^'^~^^^-^~^)'^'°^" < 1/2, 



for all n > no, since we choose A > 2 and, by assumption, nG(b) > 51ogn. Therefore, 



F[Nb >n]> 
Then, by taking the logarithm. 



4n 



(1 _ (1 _ ,)G{b)) 



n+l 



log F[Nb > n] > - log 4 + (a(l + e) - 1) log G{b) + (n + 1) log(l - (1 - e)G(6)) - log n 

> — aelogn — (q(1 + e) — 1) logn + (n + 1) log(l — (1 — e)G{b)) — logn, 

where in the last inequality we used the assumption logG(6) > log((^logn) — logn > — logn, 
for n > e*^/^, and that log 4 < aelogn for large n. Hence, 

log¥[Nb > n] > - q(1 + 2e) log n + (n + 1) log(l - (1 - e)G{b)) 

> - a(l + 2e) log n + n(l + 2e) log(l - G{b)) 
= - (1 + 2e) [a log n - n log(l - G{b))] . 

Finally, dividing by (alogn — nlog(l — G{b))) > and replacing e with e/2 yields 

logF[jV, >n] ^ ^ 

a log n — n log(l — ~ 

Symmetric arguments hold for the case where a < 1. We omit the details. 

The remainder of the proof for nG{b) < 5 logn, n > no, is deferred to Section HI □ 



2.1 Exact Asymptotics 

In this section, we derive the exact approximation for ¥[Nb > n] that works uniformly for all 
n,b sufficiently large (Theorems 12.31 and 12. 4p . Interestingly enough, this characterization is 
explicit in that it is a product of a power law and the Gamma distribution 

¥[Nb >n]^ — — — / e-'z'^-^dz, (2.6) 

^ ^ n"£(nAP[^>6]-i)y_„logP[^<,] ^ ' 

where x Ay = min(x,y) and is slowly varying. Implicitly, the argument of i{x) is altered 
depending on whether nP[A > b] < C or n¥[A > b] > C for some constant G. Hence, we 
can choose G = 1 since £(n A 1/F[A > b]) « £{n A G/F[A > b]) for large n,b. Note that 
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when — relogP[A < &] i 0, the power law dominates, whereas when — nlogP[A < 6] — )• oo, the 
integral determines the tail with the geometric (exponential) leading term. 

We would like to point out that approximation ()2.6p actually works well when F[A > b]~^ 
is large rather than b; this can be concluded by examining the proofs of the theorems in this 
section. Hence, formula (|2.6p can be accurate for relatively small values of b provided that A 
is light-tailed. This may be the reason why we obtain such accurate results in our simulation 
examples in Section [3] for b < 10. 

First, in Theorem 1 2. 3 1 we precisely describe the region where the distribution of A'^;, exhibits 
the power law behavior, nP[yl > b] < C, for any fixed constant C. Then, Theorem [231 covers the 
remaining region, n¥[A > b] > C, where PlNj, > n] approaches the geometric tail. Additional 
discussion of the results and the treatment of some special cases are presented at the end of 
this section; see Propositions 12.2 1 and 12.31 

Theorem 2.3 LetF[L > x]'^ = £{F[A > x]-^)F[A > x]'" , a > 0, x > 0, and C > be a fixed 
constant. Then, for any e > 0, there exists tiq such that for all n > uq, and nP[yl > b] <C , 

P[iVfc > n]n"£(n) ^ 
ar(-nlogP[A< 6],a) " ' ^ ' ' 

where T{x,a) is the incomplete Gamma function defined as e~^ z"~^dz. 
Proof: See Section HI 



Remark 4 Related result was derived in Theorem 3 of [19j where it was required that n < 
G{b)~^,0 < rj < 1. Note that here we broaden the region where the result holds by re- 
quiring n < C/G{b), which is larger than n < G{b)^''^. Furthermore, this is the largest 
region where the exact power law asymptotics 0(n^"/^(n)) holds since for nG(b) > C, 
T{nG{b),a) < r(C, a) ^ as C ^ oo. 



Remark 5 Note here that the incomplete Gamma function T{a,x) = z°'~^e~^dz can be 
easily computed using the well known asymptotic approximation (see Sections 6.5.32 in [IJ), 
as X — >■ oo, 

a 



1 + 



1 (a-l)(a-2) 

I o ~r . . . 



X 



X^ 



The following corollary is an immediate consequence of Theorem 
small generalization of Theorem 2.1 in [9]. 



and it represents a 



Corollary 2.1 IfF[L > x]"^ = i{F[A > x]-'^)F[A > x]~", x>0,a>0, where i{x) is slowly 
varying, then, as n ^ oo and nF[A > 6] — )■ 0, 



[Nb > n] 



r(a + i) 



(2.i 



Now, we characterize the remaining region where nF[A > b] > C. Informally speaking, 
this is the region where F[Ni, > n] has a lighter tail converging to the exponential when 
n » G{b)~^. In the following theorem, we need more restrictive assumptions for i{x); see 
the discussion before Theorem 12.21 In particular, we assume that i{x) is slowly varying and 
eventually differentiable with i'{x)x/i{x) — )• as x — )• oo. 
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Theorem 2.4 Assume that F[L > x]-^ = i(P[A > x]-^)F[A > 2;]-",a > 0, x > 0, where i{x) 
is slowly varying and eventually dijferentiable with £'{x)x/£{x) — t- as x ^ oo. Then, for any 
e > 0, there exist bQ,nQ, such that for all n > ng, 6 > bo, n¥[A > b] > C, 



F[Nb > n]n"i{F[A > b]-^) 



ar(-nlogP[A <b],a) 



1 



< e. (2.C 



Remark 6 Observe that Theorems 12.31 and 12.41 cover the entire distribution ¥[Nf) > n] for ah 
large n and b. Interestingly, the formula for the approximation is the same except for the argu- 
ment of the slowly varying part, which equals to n and F[A > b]~^, respectively. Furthermore, 
when nF[A > b] = C the formulas are asymptotically identical as £{n) = £{CF[A > b]^^) ~ 
£{F[A > as n ^ oo. 

Remark 7 Note that most well known examples of slowly varying functions satisfy the condi- 
tion — )• as X — )• oo, e.g., log^ x, log^(log x), /3 > 0, exp(logx/loglogx),exp(log'^x), 
for < 7 < 1 [see Section 1.3.3 on p. 16 in [5]]. 

Proof: Recall that 

F[Nb >n]= E[l - G(Lfe)]" 



Now, from the fact that £{x) is eventually differentiable {x > xq) and slowly varying with 
i'{x)x/£{x) ^ as X ^ oo, it follows that dF{x) = (1 + o{l))aG{x)°'-'^£-'^{l/G{x))dG{x) as 
X oo. Thus, for any < e < 1, we can select xq large enough such that 

nN. > m < (1 - CN))" - (1 + .)v^ £ (1 - cw)" °^^^;;;^^' 

/•G{xo) a—lj 

which follows from the absolute continuity of G(x), i.e., G{A) is uniformly distributed in [0,1]. 

Now, for a > 1, we consider two different cases: (a) nG{b) > logn and (b) G < nG{b) < 
log n. 

Case (a): nG(b) > log re. Observe that, for any fixed H > a + 6, we can make HG(b) 
small enough by picking bo large. Now, by continuity of G{x), there exists xq such that 
G(xq) = HG{b); we can choose xq larger than in (12. lip by picking bo large enough. Next, 
using the elementary inequality 1 — x < e~^,x > 0, we upper bound the preceding expression 

by 



F[^. > re] < e-"^(-) + ^^^^ (1 - .r^^-''' 



m 7(5(6) ' ' H^/z) 



Kl/G{b))F{b) G(6)<?<W) ^(lA) ^ ^ 

^ i{l/G{b))F{b) la^,) 
= Io + Iu (2.12) 
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where, for the third inequahty, by the slowly varying property of l{x), G{b) ^ can be chosen 
large enough such that sup(^HG{b)y^<y<G{byi KG{b)~^)/i{y) < (1 + e)^/^- 

Now, we derive a lower bound for Ii in (I2.12p . Using the monotonicity of z°~^,a > 1 and 
since < 1, we obtain 

1 _ fHG{b) 
> ^77T^G{bY~^ /_ (1 - Z^dz 



Gib)- ^ joib) 



F£(6))Y 



n + 1 ' ' ' ' I V l-G{b) J 

where in the second inequality, we use the property of slowly varying functions i{x) < for 
X large enough. Now, observe that for all x > small enough, 1 — x > e~^^, yielding 

4n 

where the last inequality follows from the fact that n/{n+ 1) > 1/2 for n> 1 and G{b)'^^^-^^ > 
Gib)"' since e < 1. We also note that 

— 71-|-1 

where we use our assumption nG{b) > C and choose H large enough. Finally, we obtain 

h > -^G(6)°e-^"'5(''). (2.13) 
An 

Now, we proceed with proving that Iq/Ii in (12.12P is negligible as n — )• oo. To this end, 
observe that 

J ^-HnGib)^ 

< 4- 



h - G(6)"e-4"G(fe) ' 

g-(H-4)nG(6)^a+l 



< 4- 



< 4e-{^^-4)nG(fe)^a+l 

< 4e-("+2) logri^a+i^ 

where we use our assumption that nG{b) > logn > 1 for n > 2, whereas for the last inequality, 
we also use the fact that H > a + Q. Thus, the preceding expression is upper bounded by 

'-^<-<e, (2.14) 

for all n > 4/e. 

Now, we upper bound Ii in (j2.12p by changing the variables z = 1 — e""/"", 

du 



F{b)i{i/G{b)) y_„iog(i-G(fc)) 

a(l + e) /■°° ^-uM_ ^-u/nY-i 

<^ =- / diXi 

- F{b)e{l/G{b)) y-„iog(i_G(fe)) n 



12 



where for the inequaUty we use e '^1'^ < 1 and extend the integral to infinity. Thus, for a > 1, 
from the preceding expression using the inequahty 1 — < x, for x > 0, we obtain 



F{b)ni{l/G{b)) J -niog{i-G{b)) 



F((i)n"('(l/G(6)) 7-„,„g(i_G,i,) 

- F(i)°''a/c(.)) ^'-"'°^<'-°'''»-">- 

Combining ()2.14p and (|2.15p . we obtain for (|2.12p that for ah n > no, 6 > 6o, 

which completes the proof after replacing (1 + e) with (1 + e)^/^. 

The remaining case (b): C < nG{b) < logn is postponed to Section HI □ 



Remark 8 Prom the preceding two theorems we observe that F[Nb > n] behaves as a true 
power law of index a when — nlogP[j4 < 6] — )• c, < c < oo, and has an exponential tail 
(geometric) when nP[A > 6] — t- oo (n >> F[A > 6]"^). More specifically: 

(i) If -nlogP[^ <b]^c, then by Theorem \T3\ as n oo, nF[A > 6] c, 

(ii) IfnP[yl > 6] — ^ oo, then — nlogP[^ < 6] — ^ oo and thus, asn — ^ oo,6 — ^ oo,nP[y4 > 6] — > oo, 

which follows from Theorem 12.41 and the asymptotic expansion of the Gamma function (see 
Remark \5\ of Theorem I2.3p . 

Interestingly, one can compute the distribution of P[A';, > n] exactly for the special case 
when the parameter a takes integer values and i{x) = 1. 

Proposition 2.2 If¥[L > x] = ¥[A > x]", for all x > and a is a positive integer, then 

F[L < b\ [a — i)\[n + z)\ 
Proof: It follows directly from (j2.10p using integration by parts. □ 

Finally, in the following proposition, we describe the tail of P[A''fc > n] for fixed and possibly 
small b. This complements the conclusion of Remark [8^ii) , however, we need i{x) = 1. 

Proposition 2.3 Let b be fixed. IfF[L > x] = F[A > x]°, a>0, x>0, then 

F Nb >n] ^ ^- ^ ~ ^ as n ^ oo. 

^ ^ F[L<b] n + 1 

Proof: See Section 4 in |8|. □ 
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10° lO' 10" 



Retransmissions: n 

Figure 2: Example 1(a). Exact asymptotics for 
a> 1. 




10° lO' 10" 10' 



Retransmissions: n 

Figure 3: Example 1(h). Exact asymptotics for 
a <l. 



3 Simulation Experiments 

In this section, we illustrate the validity of our theoretical results with simulation experi- 
ments. In all of the experiments, we observed that our uniform exact asymptotics is literally 
indistinguishable from the simulation. In the following examples, we present the simulation 
experiments resulting from 10* (or more) independent samples of A^^^^j, 1 < i < 10*. This 
number of samples was needed to ensure at least 100 independent occurrences in the lightest 
end of the tail that is presented in the figures {Ni^^i > nmax); thus providing a good confidence 
interval. 

Example 1. This example illustrates the uniform exact asymptotics presented in Theo- 
rems [2]3] and [231 i.e., approximation (12. Gh . which combines the results from both theorems. 
We assume that L and A follow exponential distributions with parameters A = 2 and /i = 1, 
respectively. It is thus clear that Fix) = e~'^^ = G{x)", where a = 2 and £{x) = 1. Now, 
approximation ()2.6p states that F[Nb > n] is given by (1 — e~^^)~^2n~^r(ne~^, 2). Note that 
we added a factor F[L < b]~^ = (1 — e~'^^)~^, as in Propositions 12.2] and [2T3l for increased 
precision when b is small; we add such a factor to approximation ()2.6p in other examples as 
well. We simulate different scenarios when the data sizes Lf, are upper bounded by b equal to 
1, 2 and 4. The simulation results are plotted on log- log scale in Fig. [2l 

From Fig. O we observe that the numerical asymptote approximates the simulation exactly 
for all different scenarios, even for very small values of n (large probabilities). We further 
validate our approximation by considering scenarios where L, A are exponentially distributed 
but a < 1; in fact, this case tends to induce longer delays due to larger average data size 
compared to the channel availability periods. In this case, we obtain a = 0.5 by assuming 
A = 1 and ^ = 2. Again, the simulation results and the asymptotic formulas are basically 
indistinguishable for all n, as illustrated in Fig. [3l 

For both cases, we deduce that for h small the power law asymptotics covers a smaller region 
of the distribution of A'^^ and, as n increases, the exponential tail becomes more evident and 
eventually dominates. As b becomes large - recall that 6 — )• oo corresponds to the untruncated 
case where the power law phenomenon arises - the exponential tail becomes less distinguishable. 

Example 2. This example demonstrates the exact asymptotics for the exponential tail 
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Retransmissions: n Retransmissions: n 



Figure 4: Example 2. Power law versus expo- Figure 5: Example 3. Power law region in- 
nential tail asymptotics. creases for lighter tails of L, A. 

as n — )• oo and b is fixed, as in Proposition 12.31 Note that this proposition gives the exact 
asymptotic formula for the region n S> 1/G{b) and lends merit to our Theorems 12.21 and 12.41 
Informally, we could say that a point rih such that — n6log(l — G{b)) nhG{b) = alogn;, 
represents the transition from power law to the exponential tail. We assume that L, A are 
exponentially distributed with A = 2 and = 1 (as in the first case of Example 1). Roughly 
speaking, we can see from Fig. |4]that the exponential asymptote appears to fit well starting 
from Uf, ae'', i.e., rif, 6, 15, 100 for 6 = 1,2, 4, respectively. 

Example 3. This example highlights the impact of the distribution type of channel 
availability periods G{x) = ¥[A > x]. We consider some fixed 6, namely 6 = 8, and assume that 
the matching between data sizes and channel availability, as defined in Theorems 12.31 and 12. 4^ 
is determined by the parameter a = 4. We assume Weibull^ distributions for L, A with the 
same index k and HL,fJ'A respectively, such that a = {^ia/i^l)^- The simulations include three 
different cases for the aforementioned distributions: Weibull with index k = 1 (exponential) 
where /Ul = 1 and /^a = 4, Weibull (normal- like) with index /c = 2 (/ul = 1,/iA = 2) and 
Weibull with k = 1/2 {fi^ = 1, fJ-A = 16). Fig. [5] illustrates the exact asymptotics from equation 
()2.6p . shown with the lighter dashed lines; the main power law asymptote appears in the main 
body of all three distributions. We observe that heavier distributions (Weibull with k = 1/2) 
correspond to smaller regions for the power law main body of the distribution ¥[Nh > n]. On 
the other hand, the case with the lighter Gaussian like distributions for k = 2 follows almost 
entirely the power law asymptotics in the region presented in Fig. [5j This increase in the 
power law region can be inferred from our theorems, which show that the transition from the 
power law main body to the exponential tail occurs roughly at n?, ~ G{b)~^. Hence, the lighter 
the tail of the distribution of A, the larger the size of the power law region. 

Example 4. In this last example, we study the case where there is a more general 
functional relationship between the distributions of availability periods A and data sizes L, as 
Theorems 12.31 and 12.41 assume. In particular, we consider the case F{x) = G{x)°' / l{G{x)~^), 
where l{x) is slowly varying. We validate the approximation (|2.6p in this more general setting. 

^In general, a Weibull distribution with index k has a complementary cumulative distribution function 
P[X > x] = e"'^/")', where ^ is the parameter that determines the mean. 
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Retransmissions: n Retransmissions: n 

Figure 6: Example Exact asymptotics for Figure 7: Example 4(^)- The asymptotes from 

the case where L follows the Gamma distribu- Theorems \2.S\ and \2.4\ for the case where L fol- 
tion. lows the Gamma distribution. 

In particular, the availability periods A are exponentially distributed with parameter ^ 
while the data sizes L follow the Gamma distribution with parameters (A,/c); the tail of the 
Gamma distribution function is defined as X^T{k)~^ e~'^^x^~^dx = r{\x, k)/T{k) and, 
therefore, the tail distribution of L can be approximated by F{x) ~ {X^~^ /r{k)) x'^^^e^'^^ for 
large x. We can easily verify that F[x) = log G'(x)~^)G(x)", where a = A//i and 

f{x) = X^'-^Tik)-^ / e-'{z/X + x)^-^dz. 
Jo 

Hence, the slowly varying function in Theorems 12.31 and 12 . 4 1 is £(x) = l/f{fi^^ logx). From the 
preceding integral representation for /(x), it can be easily shown that i(x) ~ T{k)a^^^ log^~'^ x, 
which is indeed slowly varying, and F{x) ^ {a^~^ /T{k))log{G{x)~^)''~^G{x)'^. We take A = 
2,k = 2 and /u = 2 and run simulations for b = {2, 3, 4}. In Fig. [6l we demonstrate the results 
using the approximation (12. 6p . Interestingly, our analytic approximation works nicely even for 
small values of n and b although the conditions in our theorems require n and b to be large. 

In Fig. [71 we elaborate on the preceding example. To this end, we plot two asymptotes: (i) 
the 'Initial Asymptote' corresponding to the power law asymptote provided by Theorem l2.3l and 
(ii) the 'Tail Asymptote' from Theorem 12. 4[ Combining the two, we derive the approximation 
(|2.6p . as we have already shown in Fig. [6l Hereby, we see from Fig. [7|that both asymptotes 
are needed to approximate the entire distribution well, i.e., the 'Initial Asymptote' fits well 
the first part of the distribution, whereas the 'Tail Asymptote' is inaccurate in the beginning 
but works well for the tail. Recall that these two asymptotes differ only in the argument of 
the slowly varying function £{■), which is equal to n for the 'Initial Asymptote' and G{b)~^ for 
the tail. 

4 Proofs 

In this section, we present the proofs of Theorems 12.11 and 12.31 as well as the remaining parts 
of the proofs for Theorems 12.21 and 12. 4[ 
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Proof: [of Theorem 12. 1| Without loss of generahty, we assume e < 1. We first observe that 
for n > no the assumption F[A > b] < implies that b> bo for some bo; bg can be made 

arbitrarily large by increasing uq. Hence, the upper bound follows from Proposition 12.11 
For the lower bound, we have 

F[Nb >n]= E[l - G{Lt)r 

> E [(1 - GiLbTHLb > xo)] 

> {i-G{xo)rnb>L>xo), 

since P(L < b) < 1. Now, by our assumption G{b) < l/n^'^'^ and the continuity of G(x), there 
exists xq < b such that G{xq) = 2/n. Hence, 

F[Nb >n]>(l-^) (F(xo) - Fib)) 



where we use our main assumption G(x)"(^^'^/^) < F(x) < G{x)"^^ ^^^'^ and the inequality 
1 — X > e~^^ for small x. Now, by the assumption that G{b) < 1/n^^'^, we have 



F[Nb >n]>e 



2xa(l+e/3) / ^ xa(l-e/3) 



Next, observe that (1 + e)(l — e/3) > (1 + e/3), since e < 1, and thus, 

/ 2" 1 

F[Nb >n]> e-' 



„a(l+e/3) „a(l-e/3) 

2° - 1 K 



^a(l+e/3) ^a(l+e/3) ' 

where we set = e~^(2" — 1) and by taking the logarithm, we obtain 

logP[iVfe > n]> log/ie - a(l + e/3) log n. 

Finally, since logn is increasing in n, we can choose uq such that for all n > uq, 2aelogn > 
—3 log he, i.e., 

log P[iV6 > n]> -a{l + e) log n. 

□ 

Next, we prove Theorem 12.21 in the region where n = o{logn/G{b)). 
Proof: [of Theorem 12. 2 j Here, we prove the result for nG{b) < (5 logn, for n > no, 5 > 0. 
Note that the complementary part of the proof assumes any fixed 5 > 0, and thus we can set 
5 = Sf: = ea/4:. This assumption implies that b is large, say 6 > 5o, when n is large, and thus the 
upper bound follows from Proposition 12.11 For the lower bound, it is sufficient to prove ()2.4p 
instead of ()2.5p since for b large enough, a logn < — nlog(l — G{b)) + a logn < a(l + e) logn. 
Without loss of generality, we may assume that e < 1/3. Therefore, 

F[Nb >n]= E[l - G{Lb)r 

> E [(1 - G(Lfe))n(Lf, > xo)] 

> (l-G(xo)rP[6>L>xo], 
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and since G{b) < (5 log n/n, by continuity of G{x), there exist xo,xi{xq < xi < b), such that 
G{xo) = 2(5 log n/n and G{xi) = (5 log n/n. Hence, 



[iV6>n]> (l-^^) [F{xo) - Fix,)] 

> e-^-^i^g" [G{xoT£{G{xo)~'^) - Gixi)''£{G{xi)-^)] 



where we use our main assumption and the inequality l—x>e for small x. Next, we observe 
that for all xq < x < xi, (1 — e)^(nlog~^n) < £{i6nlog'^ n) < (1 + e)^(nlog~^ n), i = 1,2. 
Thus, 



>n]>n 



-4(5 



N«/ 1 -1 ^/251ogn\" , N/,/ , 1 . f Slogn 
(1 -e)^(nlog ^n) ( - (1 + e)£(n log"^ n) ( ^ 



n 



> ^-.siin/logn)6'^log-n ^^^^^ _ _ ^ ^^^^ 

Next, since 6 = ea/4, we have 

,™r,T 1 „f£(n/logn)(5"log°n, 
F[Nh >n]> n'"'^ — — he, 

where we set he = (5°(2°(1 — e)" — (1 + e)"^) > since e < 1/3. Now, recalling the slowly 
varying property of i{n) that i{n) > n"^", for large n, we have 



a(l+2e) - „a(l+2e) ' 



since logn > 1 for all n > 2. 

And by taking the logarithm, we obtain 

log F[Nb > n]>loghe- a{l + 2e) log n 
> — a(l + 3e) logn, 

since we can choose uq such that for all n > no, aelogn > —log he- Finally, dividing by 
a logn > 0, n > 1, yields 

l2£3^>-(i + 3.), 

a logn 

which completes the proof after replacing e with e/3. □ 

Next, we present the proof of Theorem 12.31 This proof uses some of the ideas from the 
proof of Theorem 2.1 in [9j. However, it is much more involved since one has to incorporate 
the assumption nF[A > b] < G, which ensures the power law body. 

Proof: [of Theorem 12. 3| Let <I>(x) = £(x)x". Then, ^(x) is regularly varying with index a 
and, thus, for any c > 0, 

hm ' ' = c" < oo, 

and, in particular, we can choose c = e, which implies that there exists n^ such that for 
n/e'^ > Ue, 

-^<e^(-^i). (4.1) 
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Without loss of generality, we may assume that $(.) is eventually absolutely continuous, strictly 
monotone and locally bounded for a; > since we can always find an absolutely continuous 
and strictly monotone function 

= (4.2) 
[O, < x < 1, 

which for x large enough satisfies 

This implies that, for any < e < 1 and x > xq, we have 

1/$^ ((1 + e)F(x)"i) < Gix) < 1/$^ ((1 - e)F(x)-^) , (4.3) 

where denotes the inverse function of <&*(•); note that the monotonicity of <I>*(x), for all 

X > 1, guarantees that its inverse exists. To simplify the notation in this proof, we shall use <!>(•) 
to denote $*(•). Furthermore, '&^(-) is regularly varying with index 1/a (see Theorem 1.5.12 
in [5]), implying that 

1 + 



cl>-((l + e)x)~ '^^{{l-e)x), 

as X —)■ oo. Therefore, for rj^ = ri{e) = [(1 + e)/(l — e)]^^° and x large, 

V7^G{x) < l/<^>^ ((1 + e)F{x)-^) < 1/$^ ((1 - e)F(x)"^) < r/,G'(x). (4.4) 

First, notice that the number of retransmissions is geometrically distributed given the 
packet size Li,, 

FlNb >n]= E[l - G{Lb)r 

= E[l - GiLbTliLb < xo) + E[l - GiLbTMLb > xq). (4.5) 



We begin with the lower bound. For H > G and xq as in ()4.5p . we choose Xn > xq such 
that <^'*~((1 — e)F{xn)~^) = n/H, for n large, and thus, 

P[iVfe >n]= E[l - G{Lb)r 

> E [(1 - G{Lb)ri{Lb > Xn)] 



> E 



1 Y dF{x 



.X^ 1>-((l-6)F(x)-i); F(6) ' 

where we use our main assumption (j4.3p . Now, since F{b) < 1 and using the continuity of 
F{x) and change of variables z = n/^^{{l — e)F{x)~^), we obtain, 

A/<i>^({i-e)F{b)-i) ^ nJ ^>^(n/z) 



n' ^^(n/z) z"^ 



where we use that rj^Gih) > 1/$^((1 — e)F{h) ^) from (j4.4p . which holds for h large, or 
equivalently n large by our assumption nG{b) < G. Now, we consider two distinct cases: 
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If r]^nG{b) < h, where /i > is a small constant, then 



P[iVfe > n] > (1 - e) / (l - - 
Jh ^ n 

> (1 - ef^ 



^'^{n/z) z' 



dz 



nJ 



where we use the properties of regularly varying functions that for all /i < z < if and large n, 

<I>(n) 



$(ri/z) 



> (l-e)i/2^", 



and 



$'(n/z)/$(n/z) 



az 



n 



for n > H [see (|4.2p ]. Next, using 1 — x > e for 6 > and < x < xs, for n large 

enough (n > H/xs) we obtain 



¥[Nb > n] > (1 - e)3/2 
> (1 - e)3/2^ 



a 



-5H " 



$(n) 



and by choosing 6 > l/H so that e > (1 — e)^*^^, we have 



P[iVfc > n] > (1 - e) 



>(l-e) 



>(l-e) 



>(l-e)' 



>(l-e) 



$(n) 



a 



<I>(n) 



$(n) 



a 



$(n) 



nG(b) 

nG(f)) 

oo 

nG(b) 

r 

nG{b) 



h 

nG{b) 
h 



e ^z^ ^dz 



e-'z'^-^dz 



e-^z'^-'^dz 



e-'z'^^Uz 







e^^z^'-^dz 



H 



a 



e-'z'^-^dz - 2e-"H'^-^ - — 

a 



a 



{l-eY—-T{nG{b),a) 1 



where the second to last inequality follows from the approximation for the incomplete gamma 
function for large H [see Remark [5] of Theorem 12. 3| and the last inequality uses the assumption 
nG{h) < C. Now, picking H,h such that 2e-^H°'-^ + h°'/a < er{C,a), yields 

P[iV,>n]>(l-3e)^^^r(nG(6),a). 



If /i < r/,nG(6) < C, then 



P[iVfe > n] > (1 - e) 



>(l-e) 



H 



^'(n/z) n , 
e "-TTTT — r^—rrdz 



v.nG{b) ^^{n/z) z 
H 



a 



e-'z'^-^dz, 
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which follows from the regularly varying properties in the region /i/r/^ < nG{h) < z < H. For 
the preceding integral, similarly as before, we have 



H 

V.nG{b) 



> 



> 



> 



H 

nG{b) 
H 

nG(b) 
H 

nG{b) 
H 

nGib) 



e^'z^'-^dz 



e-'z'^-^dz 



r,,nG{b) 

nG(b) 
ri,nG{b) 

nG{b) 



e-'z'^-^dz 



z'^'^dz 



e^'z'^-'dz - {nGib))' 



a 



e-'z'^-^dz 



oo n° — 1 

e-'z'^-^dz - C^^ 

H a 

a 



e-^z'^-^dz-2e-''H 



nG{b) 



a 



T{nG{b),a) 1 



T{C, a) 

where we use the approximation for the incomplete gamma function for large H, that ry" — 1 — > 
4e as e ^ and nG{b) < C. Now, letting H be such that 2e-^ H""-^ + 4eC"/a < ^feT{G,a) 
yields 



[Nb>n] > 



a 



F{b)<^>{n 

Finally, since G{b) < — log(l — G{b)), we obtain 



-T{nG{b),a) 



¥[Nh > n] > (1 - e)2(l - ^) 



a 



F{b)^{n) 



r(-nlog(l-G'(6)),a), 



which proves the lower bound after replacing (1 — e)^(l — -y/e) with 1 — e. 
Next, we derive the upper bound. Note that for xq as in (|4.3p . 



(4.6) 



F[Ni, >n]= E[l - GiLb)r 

< (1 - G'(xo))" + E[l - G(L,)]"1(L, > xo) 

< e-"G'(^o) + E ( 1 



$^((l + e)F(Lf,)-i) 



- l(Lf,>xo) 



(4.7) 



which follows from (j4.3p and the elementary inequality 1 — x < e ^. Now, for any H > 
max (C, 1), we obtain 



F[Nb >n]< e-"'5(^«) + E ( 1 
< e-"^(^«) + E f 1 



1 



$-((i+.)f(l;F)) 



H 
< — 



$^((l + e)F(Lb)-i)y \^^{{l + e)F{Lf,)-^) n 



riog(n/nE)] 
fc=[log H\ 

h + h + h + h- 



n 



«^^((l + e)F(Lb) 



-1^ 



< e 



fe+i 



+ e 



-n/rie 



{4.i 
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First, we upper bound Ii in (j4.8p . which equals 

rb 



1 + e 



fl "^^^ ^^^^ 4>'(n/z) n 
Jn/'S>^((l+e)F(b)-^) ^ nJ $(n/z) $(n/z) z2 



<^{n)F{b) Jn/<s>^{{i+e)F{b)-^) ^ riJ $(n/z) $(n/z) z' 
1+e /■■'^ / <I)(n) (^'{n/z)n 

-W)mJnGib)/,y~n) ^{n/z) <^{n/z) ~^ ^' 

where we use the change of variables z = n/<I>^((l + e)F{x)^'^) and the absolute continuity of 
F{x). For the last inequality, observe that 1/$^((1 + e)F(6)"^) > G{b)/r]^ from ()4.4p . Now, 
similarly as before, we consider two cases: 

If nG{h) < rj^he < C, where /i^ > is a small constant, Ii is upper bounded by 

i,< r (i-'-Y '^^^^ zLri^ + p f ^ <h\ (4 9) 

^ - F(6)$(n) V n; $(n/z) $(n/z) z"^ + ^)F{L^,)-^) V ' 

Now, since <!>(.) is absolutely continuous and regularly varying, it follows that for all /i^ < 
z<H, 

for large re, and, by ()4.2p . 

^'(n/z) 
$(re/z) re ' 

for n > H. 

Next, we compute the second term in (14. 9p as 

F(L)< ' + ' 



^ 1 + g ^ (1 + 6)2/1^ 



$(re//i,) 



F(6)$(re//i,) - $(re) ' 

which follows from the uniform distribution of F{L) and using $(n)/$(re//i,) < (l + e)i/2/ja 
for large re, along with F{b)^^ < (1 + e)^^^. Now, observe that the first term in ()4.9p is upper 
bounded by 



^C'^) A, ^ ^("') JnGib)/v. ^ n 

since nG{b) < h^rj^. Also, for the integral we obtain 

JnG(b)/r,^ ^ IT"' JnG{b) ^ nJ JnG{b)/r,, ^ riJ 



< I [I - -]" z'^-Uz + I " z'^-^dz 
JnG(b) ^ riJ JnG(b)/rj, 

< r fl--Vz"-idz + (reG(6))°(l-777")/a 

JnG{b) ^ riJ 

^ r. (l- -Y z^-^dz + bCe/a, 



nG(b) 



re 
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after observing that 1 — ry^ " — )• 4e as e — )• 0. Now, by changing the variables 1 — zjn = e 
we have 



h < 



< 



< 



< 



< 



a{l + e 
$(n) 



z 



(l + e)25C7"e ^ (1 + e)2/i^ 



<I>(n) 



<I>(n) 



2 /--H" , 

(1 --Y z'^-^dz 
nG{b) ^ 

nlog(l-G{6)) 



$(n) 



nlog(l-G{b)) 



2 />oo 



a(l + e) 
<I>(n) 

a(l + e)^(l + yi) 
<I>(n) 



-nlog(l-G(b)) 



1 + 



$(n) 
5C°e + 



Q!r(2C,a) 



(4.10) 



»^log(l-G(fe)) 



where, in the second inequahty, we use e~"/" < 1, the inequahty 1 — < x, x > and 
extend the integral to infinity. Last, we pick e small, such that 5C"e + hf < y/€aT{2C,a). 

Note that the preceding equation along with (]4.6p imply that Ii > {l—€)aT{2nG{b), a)/^{n) 
> {a/2)T(2C,a)/^{n), for ah n > uq and e < 1/2, by the inequality 1 — x > e~^^ for x > 
small, since by assumption nG{b) < C, i.e., G(b) is small. 

If her]e < nG{h) < G, we have 



h < 



1 + e 



F{b)<^>{n) 



H 



nG(b)/-n, 



1 



z\" <i>(n) <I>'(n/z) n 



nJ ^{n/z) ^{n/z) z' 



dz. 



and, by the properties of regularly varying functions in the interval n/H < n/z < 1/G(b) < 
n/he, for H > G, and using the same arguments as in (j4.10p . we have 



,2 rH 



h < 



< 



nG{b)/r,, 



< 



a{l + e 
<l>(n) 

a(l + e) 
<I>(n) 

a(l + e)2(l + Vi) 



2 fOD 



-) z^-^dz 

71/ 

e^^z'^^^dz 



-nlog(l-G(fe)) 



1 + 



ar(2C, a) 



f 



e-'z^-^dz. 



^("') J~~nlog(l-G(b)) 

Therefore, from both cases, it follows that for all n > hq 

a{l + e 



h < 



$(n) 



-r(-nlog(l-G(6)),a), 



after replacing (l + e)^(l + ^/e) with 1 + e. 

Next, we evaluate the second term in (j4.8p as 



|"log(7l/ne)] 

l2= E e- 

fc=Liog Hj 

riog("/"e)i 

fc=[logHJ 



n 



ci>-((i + 6)F(L,)-i; 



n 



< e 



n 



(4.11) 
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< 

k=llogH\ 
oo 

_ fc i + e 



F(6)$(n/e'=)' 

fe=[log/Jj 

which follows from the fact that the integral in the second inequality is equal to ¥[F{L) < 
(1 + e)/^> (n/e'=)]/F(5) and F{L) is uniform in [0, 1]. Thus, 

^ - F(b)^(n) ^ ^inle^) 
- F(5)^>(n) ^ 

where we make use of (|4.1|) . Since the preceding sum is finite, we obtain that for large H and 
all n > no, 

/2 < (4.12) 
Last, we observe that, for fixed xq, it follows that for n > uq, 

Jo + /3 = e-^'^C^o) + g-n/n, < 1/^. (4.13) 

Finally, using (I4.11l) - (j4.13p . we obtain for (|4.8p that for all n > no, 
mb>n] < (l + e)2_^r(-nlog(l-G(6)),a), 

which completes the proof after replacing (1 + e) with (1 + e)^/^. □ 

Proof: [of Theorem I2.4j Here, we complete the proof of Theorem 12. 4t recall that a > 1. 
Case (b): C < nG{b) < logn. In this region, for any fixed H > 5, we choose the smallest 
m > 1 such that iJe™ - 4 > (a + 2) log n/C, i.e., m = [log ((a + 2) log n/C + 4) - log H~\ . 
Furthermore, it is important to note that this choice of m allows for He"^G(b) to be small 
enough, since He"^G{b) < He"^ log n/n = 0(log^ n/n) — )• 0, as n — )■ oo, by the assumption that 
nG{h) < logn. Then, by continuity of G{x), there exists xq such that G{xq) = He^G{h) (larger 
than Xq in (|2.1ip for 6o large enough) and using the elementary inequality 1 — x < e~^,x > 0, 
we upper bound the expression in (12. IIP by 



^[N, > n] < e-"^(-') + + 

!< [b) 

F{b) 



k=0 



Gib) ^ ni/z) + ^„ JHe>^Gib) ^ 



:i - zTz'^-^dz 



^(1/G(6)) Jcub) 

m-l ^ Nl/2 _ [■He''+^G(b) 



He^G{b) 



where the last inequality follows from the monotonicity of ^ for a > 1 and the slowly varying 
property of l{x), svLY>(^HG{b))-'^<z<i/G{b)) / (b)) / i{z) < (1 + e)^/^, while for the second term. 
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we have that sui)(^He'=+^G{b))-^<z<{He''G{b))-^ K'^/ie''G{b))/i{z) < (1 + e)V2^ k = . . . {m - 1), 
since He"^G{b) is smah enough. Now, since i{x)/£{x/e) < e for x large, it follows that 



F(6)£(l/G(6)) 



^ F{b)ni{l/G{h)) ^ ^ ^ 



k=0 

A 



h + h + h- (4.14) 

Now, we derive a lower bound for Ii following similar arguments as in (I2.13p . Note that, for 
a; > small enough, 1 — x > e~^^, and thus, for H large enough, we have 



m/cwm,. > e,,).-. r"V 

a(l + e) JG{b) 

> (1 - ^(6))"+! - (1 - HGlb)r+^ 



n+l 

\n+l 



_™«-i(i^G(6)r^ I ( i-HGib) y^'\ 

-""^'^ n + l ['-[-T^gW) ) 

> 9E^e'~^nGi'\ (4.15) 

where the expression in brackets is bounded from below by 1/2 as in (j2.13p . 

Now, we prove that Iq/Ii in (j4.14p is negligible as n — )• oo. To this end, observe that 

To ^ 4F(6) e-^^'""^W^(l/G(fe))n 
h - a(l + e) G(6)"-ie-4"G{6) ' 

where we use (|4.15p . Next, since a > l,F{b) < 1, and using the slowly varying property that 
i{x) < X for large x, we obtain 

h ~ G{b)'' 
and since nG{b) > C, we have 

lo e-(J^e'"-4)nG(fe)^a+l 

A (nG(6))" 

< 4C-°e-(^^'"-^)^n-+i 

< 4C-"e-("+2)logn^a+l^ 

where the last inequality follows from the fact that m was chosen so that (He"^ — 4) > (a + 
2) log n/C. Thus, the preceding expression can be rewritten as 

In 4 

< ^ < e/2, (4.16) 

for all n > SC^^/e. 
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Next, for the ratio I2/I1 we proceed similarly as before 



m— 1 



< 4 ^ gfcg-(He'=-4)nG(fe)(^_g-gfe+l^a-l^ 

fe=0 

m— 1 

< 4i?""^ ^ gfcg-(He'=-4)Cga(fc+l)-fc-l 

fc=0 

00 

< AH'^-^e'^^ Yl e-5(^'-i)^+°(*^+i)-^ < e/2, (4.17) 

fe=0 

where for the last inequality we use H > 5. Now, we further observe that the preceding sum 
is finite and thus, letting H — )• 00, the above ratio converges to 0, i.e., I2 < e/i/2 for large H. 

Hence, since the upper bound for Ii from (|2.15p holds in this case as well, by putting (j4.16p 
and (j4.17p together, we obtain for (j4.14p that for all n > no, 6 > 

which completes the proof after replacing (1 + e) with (1 + e)^/^. 

Last, we prove the lower bound for nG{h) > C; here, we do not need to distinguish two 
cases. Thus, starting from ()2.10p and proceeding with similar arguments as in the proof for 
the upper bound, we obtain 



a(l - e) r-in+l)log(l-HG{b)) -um _ ^-^^^a-i 



F{b)i{l/G(b)) y_(„+i) iog(i-G(6)) ~ n + l 

where we use the slowly varying property of i{x) and pick xq < b such that G{xo) = HG{b). 
Next, using the inequality 1 — > {1 — 6)x, for some 6 > and all x > small enough, we 
have 

Q,n _ gVl _ A^a-l l—nlog{l-HG{b)) 

^ ' J- F{b)l{l/G{b))n L„iog(i-G(6)) ^n) 

^ «(l-g)^ /"^"^("^ 
> — — — - — f-— - / e u du 



F{b)n-e{l/G{b)) y_„iog(i-G(6)) 



A 



Fib)n^i{l/G{b)) 
h-h, 



nlog(l-G(fe)) JHnGib) 



where, in the second inequality, we choose 6 > small enough such that (1 — 6)" ^ > (1 — e) 
and that —n log{l — HG(b)) > HnG{b). Next, we proceed with showing that I2/I1 is negligible 
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for large n. Note that — nlog(l — G(6)) < 2nG{b), which follows from the elementary inequality 
g-2x ^ I _ £qj. all X > small enough. Thus, 

" JZ&b) e-'^u'^-^du 



< 



2nG(b) 

(2nG'(6))"-ie-2"G(6) 



where we use the approximation for the incomplete gamma function for large H [see Remark 
of Theorem 12 .Sj . Now, using the main assumption nG(b) > C, we obtain 

^-^ < 2F-ie-(^-2)^ < e, 
h 

for H large enough. Then, using the preceding observation, we obtain 

which completes the proof after replacing e with e/3. 

Now, if a < 1, the proof uses almost identical arguments coupled with the fact that u°'~^ 
is a decreasing function. We omit the details to avoid unnecessary repetitions. □ 



A A note on the uniform distribution of F{L) 

Here, we give a brief comment on the fact that F{L) is uniformly distributed if L is a continuous 
random variable, i.e., F(x) is absolutely continuous. This statement is immediate when F{x) 
is strictly increasing, e.g.. Proposition 2.1 in Chapter 10 of [T7], since its inverse F^(x) is well 
defined and F{F^{x)) = F*~{F{x)) = x, meaning that, for < x < 1, 

¥[F{L) <x] = F[L < F^{x)] = F{F^{x)) = x. 

However, when F{x) has flat intervals, F'^{x) is not uniquely determined, e.g., we may use 
a standard generalized inverse: F^{x) = inf{y : F{y) > x}. Now, we still have F{F^(x)) = x, 
but F'^{F{x)) 7^ x, in general; the equality only holds if x is a point of increase for F{x). But, 
since F{x) is absolutely continuous, we can exclude the flat regions of F{x). Formally, the 
derivative f{x) = F'{x) exists a.e.-Lebesgue, and where it does not, we can set f{x) = 0. 
Furthermore, F{x) = f{u)du. Thus, since F'^{x) is strictly increasing, for < x < 1, 

F[F{L) <x]= F[F^{F{L)) < F^{x)] 

/oo 
l{F^iF{u)) < F^{x), f{u) > 0)dF{u) = F[L < F^{x)], 
-oo 

since F{F^{u)) = u when f{u) > 0, implying that F[F{L) < x] = F{F^{x)) = x. 
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